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ABSTRACT 

Potential tests using turbulence velocity histories measured 
in the SPAN-MAT program are outlined to determine validity of the 
homogeneous, momentarily frozen assumptions for the vertical 
turbulence velocity component and the homogeneous, isotropic, 
momentarily frozen assumptions for the horizontal turbulence 
velocity components. In addition, methods are reviewed -for pre- 
diction of the crosscorrelation function between any two 
spatially separated turbulence velocity components using the 
homogeneous, isotropic, momentarily frozen assumptions and 
measurements of the transverse and longitudinal turbulence 
velocity components. 
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INTRODUCTION 

This memorandum outlines potential tests using turbulence 
velocity histories measured in the SPAN-MAT program to determine 
validity of the homogeneous/ momentarily frozen assumptions for 
the vertical turbulence component and the homogeneous/ isotropic/ 
momentarily frozen assumptions for the horizontal turbulence 
components. In addition, methods are reviewed for prediction of 
the crosscorrelation function between any two spatially separated 
turbulence components using the homogeneous, isotropic, momen- 
tarily frozen assumptions and measurements of the longitudinal 
and transverse turbulence components made using a single probe. 

Turbulence measurements are being carried out in the SPAN- 
MAT program using a B57-B aircraft equipped with probes at three 
separate locations. There is a nose probe and a probe located 
under each wing tip on a pylon designed to accept fuel pods. The 
wing-tip probes are about 60 feet apart. All three probes lie 
approximately in the same (horizontal) plane. The nose probe is 
substantially forward of the two wing-tip probes. The probe 
configuration is illustrated in Figure 1, where the two wing-tip 
probes are labeled a and c, and the nose probe is labeled b. 

Each probe is capable of measuring three components of turbulence 
velocities: the forward component u, lateral component v, and 

vertical component w. 

The directions of the turbulence velocity components u, v, 
and w are oriented with respect to the direction and position of 
the measurement aircraft flight path. In describing turbulence 
correlation measurements, one also must be concerned with the 
orientation of the straight line connecting the positions of the 
two measurement points relative to the direction of the turbu- 
lence velocity vector of the component of turbulence being 
measured. In this memorandum, two such components are dealt 
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with: the longitudinal and transverse velocity components. The 

directions of these two components relative to the two measure- 
ment locations are illustrated in Figure 2. 
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LONGITUDINAL TURBULENCE VELOCITY CORRELATION 



TRANSVERSE TURBULENCE VELOCITY CORRELATION 


Pig. 2. Longitudinal and transverse turbulence velocity 
correlations . 
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TESTS OF HOMOGENEOUS MOMENTARILY FROZEN HYPOTHESIS 
FOR THE VERTICAL (TRANSVERSE) TURBULENCE COMPONENT 

On physical grounds » we would expect the vertical component 
of atmospheric turbulence to be statistically spatially homo- 
geneous over areas in a horizontal plane large in comparison with 
the dimensions of an aircraft. The simplest test of such spatial 
homogeneity of the vertical turbulence component is to measure 
and compare the autocorrelation functions of the vertical 
turbulence component obtained from each of the three SPAN-MAT 
probes. In reality, such a comparison probably is more of a test 
of the uniformity of the probes and associated electronics than 
it is of spatial homogeneity. In particular, the rms velocities 
of the vertical turbulence components obtained from each of the 
three probes should be compared. 

A second simple test employs both the spatially homogeneous 
assumption and the momentarily frozen assumption (Taylor's 
hypothesis). Denote by angular brackets <•••> a long time 
average ; i .e . , 

A 1 t ' +T 

<f(t)> = i / f ( t )dt , (!) 

t ' 

which for very large T is independent of both t' and T if f(t) is 

a stationary random function. Denote the vertical turbulence 

component measured at each of the three probes a, b, and c in 

Figure 1 by w a , w b , and w c , respectively. Denote the distances 

between probes a and b, b and c, and a and c by r ab , E’bc' an< ^ 

r , respectively. Denote the autocorrelation function of the 
d C 

vertical turbulence component as a function of a generic spatial 
separation r in a horizontal plane by R ww ( r )* Then, if the 
vertical turbulence component is assumed spatially homogeneous in 
a horizontal plane, we have 
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R ( r , ) = 
ww ab 


<w q ( t ) w b ( t ) > 


( 2 ) 


R (r, ) = 

ww be 


<w ( t )w ( t ) > 
b c 


and 


(3) 


R ( r ) 
ww ac 


<w (t)w (t)> . 
a c 


(4) 


Denote the aircraft speed by V. Furthermore, if the spatial 
homogeneity and momentarily frozen assumptions are both employed, 
we also have 


R ( r ) = 
ww ab 


<w ( t ) w (t + T )> 
a a ab 


ab 


= r /V 
ab 


(5) 


R 

ww 


(r. ) 

be 


<W ( t ) W (t+T )> , 

a a be 


and 


R (r ) 
ww ac 


<W (t)w (t+T )> , 

a a ac 



r 



L 

1 = 
ac 


r 

ac 


/V . 


( 6 ) 


(7) 


Each of the relations (5) through (7) describes the auto- 
correlation function of the vertical turbulence component 
measured at the same sensor "a" at a time delay t computed from 
the aircraft speed V and the distance r between two of the 
sensors. Sets of relations analogous to equations (5) through 
(7) can be written by replacing each w a by w b or w c . 

The homogeneous and momentarily frozen assumptions of the 
vertical turbulence component can be tested by comparing the 
crosscorrelation measurements described by the right-hand sides 
of equations (2), (3), and (4) with the autocorrelation 
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measurements described by the right-hand sides of eguations (5)/ 
(6), and (7), respectively. 

The above test can be generalized using the method suggested 
in Figure 3. Employing the homogeneous and momentarily frozen 
assumptions, for any value of r a , we have 

R (r ) = <w (t)w (t+x )> ( 8a ) 

ww a a a a 

= <W (t)W (t+T )> ( 8b ) 

b b a 

= <w (t)w (t+T )> , (8c) 

c c a 

where 

x = r /V . <9) 

a a 

Furthermore, for any value of r & > r Qc , we also have 

R (r ) = <w (t)w ( t + t ' ) > t 
ww a a c a 

where from Figure 2, it follows that 


r 2 = r 2 + V 2 ( t ' ) 2 
a ac a 


or 


a 



r > r 
a ac 


(ID 


If the homogeneous and momentarily frozen assumptions are satis- 
fied by the vertical turbulence component, for any value of 

r > r , the autocorrelation function computed by the right-hand 
q 3 C f 

side of equation (8a), (8b), or (8c) with T a given by equation 

(9) should equal the crosscorrelation function computed by equa- 
tion (10) with t' is given by equation (11). 

d 
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TESTS OF ISOTROPIC MOMENTARILY FROZEN HYPOTHESIS FOR THE 
HORIZONTAL (LONGITUDINAL AND TRANSVERSE) TURBULENCE COMPONENTS 


Homogeneity in the horizontal plane requires that the 
"forward" turbulence velocity component u pointing in the 
direction of aircraft motion have the same statistical properties 
at each of the three probe locations a , b, and c. In particular, 
this assumption implies that each of these three measurements 
u a (t), u b (t), and u c (t) possess the same autocorrelation 
function. Homogeneity also implies that each of the lateral 
components v Q ( t ) , v b ( t ) , and v c ( t ) also possess the same 
autocorrelation function, which generally will differ from that 
of the u's. 

Furthermore, isotropy in the horizontal plane requires 

directional independence of the character of the turbulence. In 

particular, isotropy requires that <u 2 (t)> = <v 2 (t)> with 
r 3 s 

analogous relations for the other two probes b and c. 

The configuration of the probes permits other tests of 
isotropy in the horizontal plane utilizing the momentarily frozen 
assumption. The autocorrelation function R LL<r ac ) of the 
longitudinal turbulence component with spatial lag r ac can be 
evaluated directly from the lateral turbulence velocities v a (t) 
and v c (t ) illustrated in Figure 4 by 


R LL (r ac ) 


<v (t)v (t)> . 

a G 


( 12 ) 


Furthermore, utilizing the isotropic and momentarily frozen 

assumptions, r ll( r ac } also can be evaluated from the forward 
turbulence component at any of the three probes, as illustrated 

in Figure 4 for the middle probe b; i.e.. 
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Fig. 4. Configuration for testing isotropic and momentarily 
frozen assumptions for the longitudinal turbulence 
component in the horizontal plane. 
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R TT (r ) = <u, (t)u, (t + t )> 
LL ac b b ac 


" < u a ( t )u a ( t + T ac )> 


= <U c (t)U c (t+T ac )> , 


(13a) 

(13b) 

(13c) 


where 


x = r /V • U 4 > 

ac ac 

Thus, the isotropic and momentarily frozen assumptions in the 
horizontal plane can be tested by comparing experimental measure- 
ments described by the right-hand side of equation (12) with 
those described by the right-hand sides of equations (13a), 

( 13b ) , or ( 13c ) . 

A similar set of relations can be written for the auto- 
correlation function of the transverse or normal turbulence 
component as illustrated in Figure 5. For the autocorrelation 

function of the transverse turbulence component RNN^ r ac^ 
evaluated at the spatial lag r ac , we have from Figure 5, 


R mnt ( r = „ ) 

NN ac 


<u (t)u(t)> , 
a c 


(15) 


which involves the forward turbulence components measured at a 
and c with no temporal lag. From Figure 5, we also see that 


R MM (r- ) 
NN ac 


<v b (t)v b (t+T ac )> 

<v a (t)v a (t+T ac )> 

<v c (t)v c (t+T ac )> , 


(16a) 

(16b) 

(16c) 
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Pig. 5. Configuration for testing isotropic and momentarily 
frozen assumptions for the transverse turbulence 
component in the horizontal plane. 
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where t is given by eguation (14). Each of the equations (16) 

3. C 

employs the momentarily frozen assumption. Isotropy and homo- 
geneity in the horizontal plane require that turbulence measure- 
ments described by the right-hand side of equation (15) equal to 
those described by each of equations (16a), (16b), and (16c). 

The assumptions of homogeneity and isotropy, and the 
momentarily frozen assumption, all in the horizontal plane, imply 
certain other relations among correlation measurements of the 
forward and lateral turbulence components. First, we note that 
the longitudinal and transverse autocorrelation functions R L L(r) 
and res P ec t ively , can be measured for any value of r from 

the forward and lateral turbulence components, u(t) and v(t), of 
any probe. For example, using probe b, we have 

R LL (r) = <u b (t)u b (t + x)> , t = r/V (17) 


and 


R NN (r) = <v b (t)v b (t+T)> , t = r/V , (18) 

where Figure 4 is helpful in understanding equation (17) and 
Figure 5 is helpful in understanding equation (18). Using Rll( r ^ 
and R NN (r), obtained as above, together with the homogeneous, 
isotropic, and momentarily frozen assumptions, crosscorrelation 
functions of the various turbulence components can be predicted 
to compare with their measured values. For example, using a 
well-known relationship [1], we have for the crosscorrelation 
functions of the forward component u and the lateral component v, 
at probe locations a and c, respectively, 
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R u u (r) ' l*LL (r) ' R NN (r) ] + %N (r) 


a c 


= 7^ R LL (r) + ^ 1 ‘ 77) R NN (r) 


r 2 r 2 

_£ R (r) + R (r ) 

r 2 LL NN 


(19) 


a c 

E v V (r) ' [R LL (r) ' R NN (r) J ? + R NN (r) 


a c 


ac 


r2 R LL (r) + Cl - %<r> 


r 2 r 2 

7 T ^ r + — 7 ^MN ^ ^ 9 

r LL r z NN 


( 20 ) 


and 


R (r ) 
u v 
a c 


r r 
c ac 


* [ R LL (r) - R NN (r) l -r* 


( 21 ) 


where the turbulence components and relevant geometry are illus 
trated in Figure 6. There is a temporal delay of the measure- 
ments u c and v c relative to the measurements u a and v a by the 
amount 


T 

C 


r 

c 

V - 


■ 2 _ 


, 1/2 


ac 


( 22 ) 


which follows directly from Figure 6. 
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Forward 



Fig. 6. Configuration for testing isotropic and momentarily 

frozen assumptions for turbulence velocity components in 
the horizontal plane. 
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Using measurements of Rll(t) an< ^ R NN^ r ^ obtained as in 

equations (17) and (18), the crosscorrelation functions R^ u (r), 

a c 

R (r), and R (r) can be predicted by equations (19) through 
v v u v 

a c a c 

(21), respectively, for any value of r > r a c’ These predictions 
can be then compared with direct measurements of these cross- 
correlation functions to test the homogeneous, istropic, 
momentarily frozen assumption, where the direct crosscorrelation 
measurements are 


R U u (r) = <u a (t)u c (t+T c }> ' 
a c 


(23) 


R (r ) - <v (t )v (t+x ) > , 
v v a c c 

a c 


(24) 


and 


R (r) = <u (t)v (t+T )> , 

u v a c c 

a c 


(25) 


where is given by equation (22)* 

An alternative to the above approach is to use equations 

(19) and (20) to solve for Rll^) and R NN^ r ^ termS of 

R (r) and R (r). The result is: 
u u v v 

a c a c 


R LL (r) 


r 2 R (r) 
c u u 
a c 


- r 2 R (r) 
ac v v 
a c 


■2 - r T 


ac 


and 


R NN (r > 


r 2 R ( r ) 
c v v 
a c 

- 2 - 


- r 2 R (r) 
ac u u 
a c 


_ ;t 


ac 


(26) 


(27) 
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provided that r c * r ac’ values of RLL(r) and f° r 

r > r predicted by ecuations (26) and (27) from the cross- 

3 C ^ 

correlation measurements (23) and (24) then can be compared with 
the values of Rll ( r ) and R NN (r) measured as in equations (17) and 
(18). However, this latter procedure might be somewhat ill- 
conditioned [ 2 ] . 


17 



Technical Memorandum No. 938 


BBN Laboratories Incorporated 


PREDICTION OF CROSSCORRELATION FUNCTION 
BETWEEN ANY PAIR OF TURBULENCE COMPONENTS 

When the homogeneous, isotropic, and momentarily frozen 
assumptions are employed, it is possible to predict the cross- 
correlation function between any pair of spatially separated 
turbulence velocity components from the longitudinal and trans- 
verse autocorrelation functions Rll( r ) and respectively, 

measured at a single probe location. Let us define, for any 
spatial location, the three turbulence velocity components by 

u x = u , u 2 = v , u 3 = w (28) 

where u is the forward component, v is the lateral component, and 

w is the vertical component. Let x denote a given location in a 

3-dimensional frozen turbulence field, and x + r denote a second 
location in the field, where r is the 3-dimensional vector 
connecting the two locations. Let r denote the length of the 
vector r ; i . e . , 

r = /rf + r| + r* , (29) 

where r x is the component of the vector r in the Uj = u = forward 
direction; r„ is the component of r in the u 2 = v = lateral 

direction; and, r 3 is the component of r in the u 3 * w = vertical 

direction. Let 

R. (r) - <u.(x) u (x+r)> ( 3 °) 

j £ ~ j~ £.~~ 

denote the crosscorrelation function between any two turbulence 
velocity components Uj (x) and u £ (x+r) at the above-mentioned two 
locations. For example, 
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R ,(r) = <u(x) v(x+r ) > 

is the crosscorrelation function of the forward turublence 
velocity component u at location x and the lateral turbulence 
velocity component v at location x+r. Similarly , 


(31) 


R~~(r) = <w(x) w(x+r)> 

3 <3 , 


(32) 




is the crosscorrelation function of the vertical turbulence 
velocity component at location x and location x+r. 

The general form [1] expressing R. (r) in terms of the 
longitudinal and transverse turbulence velocity autocorrelation 


functions is 


V r J = [R LL (r) ' R NN (r) J r 


r .r 

] — — 7 -— + R M J r > « ■ „ ' 

J NN j SL 


(33) 


where the distance r is defined by equation (29), rj and r £ ; j ,£ 
1,2, or 3 are described in the text following equation (29); and, 

6 . is Kronecker's delta; i.e., 

«. - P' 3 ' * . (34) 

11 I 0, j * t 

The geometry relevant to equation (33) is sketched in Figure 7. 

When both components j and £ in equation (33) represent the 
vertical component Uj = w, we have j = £ = 3; hence, 6 j £ = 1 . 

When the positions x and x+r of both vertical components lie in 
the horizontal plane (normal to the direction of U 3 = w), we have 


r = r = r, = 0. In this case, equation (33) reduces to 
j £ 3 

R 3 3 ( £ ) = t ) — ^ ) ' 


WW ~ 


NN’ 


( 35 ) 


which is the case discussed in equations (8) and (10). 
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When j and £ in equation (33) both represent the forward 
component Uj = u; i.e., j = £= 1/ we also have 6^ = 1, and 
r. = r = n. In this case, equation (33) reduces to 

3 £ 1 


R ll<E> 5 R UU^' = [R LL (r) - R ™ (r) ] + R NN <r) ' 


(36) 


which is the same as the first line in equation (19) when r j is 
set equal to r c « 


Similarly, when j and £ in equation (33) represent the 

lateral component U 2 = v; i.e., j = £ = 2, then 6 j ^ = 1 and 

r = r,. In this case, equation (33) becomes 
0 ^ 


R 22 (r) = 


R (r) 
vv ~ 


* t R LL (r) * 


R NN (r) ] 


R (r) 
NN 


(37) 


which is the same as the first line in equation (20) when r 2 is 
set to r ac ♦ When j = 1 and 1 = 2, we have 6 j £ = and equation 
(33) reduces to 


R 


12 <r> = R JV * t R LL (r) ' R NN (r) ) - 


r i r 2 

T~ 


(38) 


which is the same as equation (21) when r j = r c and r 2 = r ac * 

The SPAN-MAT configuration allows the longitudinal auto- 
correlation function R LL (r) to be measured at each of the three 
probes using the forward turbulence velocity component u; i.e., 


R ( r ) = R (r) = <u(t)u(t+x)> , t = r/V (39) 

Lj j_i U 

where V is the aircraft speed, and both u(t) and u(t+t) are 
measured at the same probe. The transverse autocorrelation func- 
tion r nn ( r ) can be measured using the lateral turbulence velocity 
component v(t) at any of the probes; i.e.. 
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R MM (r) = R (r) = <v(t)v(t+x)>, x = r/V (40) 

NN VV 

or using the vertical turbulence velocity component w(t) at any 
of the probes 

R. 1M (r) = R ( r ) = <w(t)w(t+x)> , x = r/V (41) 

NN ww 

where both v(t) and v(t+x) in equation (40) are measured at the 
same probe, and both w(t) and w(t+x) in equation (41) are 
measured at the same probe. 

When one or both of the two components j or £ in equation 
(33) represents either v or w, but not j representing v, and £ 
representing w, or j representing w and £ representing v, either 
the measurement (40) or (41) should be used depending on which 
component (v or w) is represented in equation (33). When j 
represents v and £ represents w, or j represents w and £ 
represents v, then judgement will have to be used as to which of 
the two measurements, equation (40) or equation (41), should be 
used to represent ec ? ua tion (33). 

To express equation (33) as a crosscorrelation function with 
temporal (rather than spatial) lag, Taylor's hypothesis must be 
used for the particular geometric configuration of interest. For 
the configuration where both x and x+r are in the horizontal 
plane on a line normal to the fuselage centerline, the trans- 
formation to temporal lag becomes 


where V is the aircraft speed and x is the temporal lag. The cross- 
spectrum of the two turbulence velocity components is obtained by 
forming the Fourier transform of equation (33) with respect 
to t after employing Taylor's hypothesis as in equation (42). 
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